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Abstract 

We reexamine the massive graviton dark matter scenario (MGCDM) which was recently consid- 
ered as an alternative to dark energy models. When introducing the native and effective equations 
of state (EoS), it is shown that there is no phantom phase in the evolution toward the far past. 
Also we show that the past accelerating phase arises from the interaction between massive graviton 
and cold dark matter. 
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1. INTRODUCTION 



w DE (a) = -1 - - = -1 - - ( ) , (1) 



Recently the massive graviton dark matter scenario (MGCDM) which has been originally 

proposed by Visser [l| was considered as an alternative dark energy model to explaining the 

present accelerating universe {2, 3|. In order to explain the evolution of the universe, they 

used the "geometric" dark energy EoS mainly as [3j 

1 d\n5H 2 2 (7- 10a 2 
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which shows that the massive graviton theory can be treated as an additional effective 
fluid with EoS ([1]) cosmologically. Also they have emphasized that the current Hubble 
function has only two free parameters of H and Q m which are the same number of pa- 
rameters as the ACDM model. However, this MGCDM showed a different phantom phase 
of lim a _j,o wde(o) — —5/2 and a current decelerating phase 3J], in compared to the ACDM 
model with = —1. This is mainly because they used the geometric dark energy EoS 
w DE (a). 

On the other hand, it is well known that if the conservation law is modified due to the 
presence of other matter as (jSj), one has to introduce the effective EoS parameter u>° ff (a) as 
well as the native EoS w° at (a) to describe the evolution of the universe correctly Q. For 
the holographic dark energy model, two of authors have clarified that although there is a 
phantom phase when using the native EoS [5fl, there is no phantom phase when using the 

n q 

effective EoS |a,LZ|. 

In this work, we wish to reexamine the evolution of the universe based on the MGCDM 
by introducing the native and effective EoS. It is hard to derive any phantom phase when 
using the effective EoS Wg S instead of the geometric EoS kjde- Finally, we show that the 
past accelerating phase arises from the interaction between massive graviton and CDM. 



2. MASSIVE GRAVITON 

We briefly review the Visser's massive gravity approach [lj. The action is given by 
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where the first term is the Einstein-Hilbert Lagrangian and the last is the Lagrangian of 
matter. The second term is designed for the massive graviton expressed in terms of the 



bi-metric (go, g) as 
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where m 5 is the graviton mass and (<7o)^ is a general flat metric. 
The Einstein equations takes the form 
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where G^ v is the Einstein tensor, T^ v is the energy-momentum tensor for a perfect fluid, 
and 



(5 - £o)<x P - ^(^o)<xp(^ ^ x (9 ~ 9o)*P 
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In the limit m g — > 0, we recover the Einstein equation is recovered. Importantly, since the 
Einstein tensor should satisfy the Bianchi identity V„(j m1/ = 0, the non- conservation law of 
the energy-momentum tensor is obtained as 

which will play the crucial role in the cosmological evolution. 



3. COSMOLOGICAL EVOLUTION 

In order to apply the action ([2]) to cosmology, we need to introduce the bi-metric explicitly. 
First we use the Friedmann-Roberston- Walker metric as the dynamical one 

d s 2 FRW = c 2 dt 2 - a 2 (t) [dr 2 + r 2 dtt 2 2 ] (7) 

with a(t) the scale factor. The flat metric expressed in terms of spherical coordinates is 
proposed to be the static one: 

ds 2 = c 2 dt 2 - [dr 2 + r 2 dQ 2 2 ] . (8) 
The evolution of a MGCDM cosmology is governed by two Friedmann equations: 
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where p is the energy density and p is the pressure. From the observation of Eqs. fl9]) and 
( FlOl) . one can read off the energy density and pressure for massive graviton as 
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We note that the positive energy density and negative pressure is allowed to explain the 
accelerating universe. In this case, it requires 



a 2 < 1. 



(13) 



Hence the native equation of state is simply given by 
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whose limit of a — > is zero. From ([H]), we find the the non-conservation law for the matter 
including cold dark matter (CDM) 
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with H = a/a the Hubble parameter and 



(16) 



On the other hand, one can obtain the non-conservation law for massive graviton by com- 
puting (flT!) and (fT2l) directly 



Pg + 3H(Pg + 
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Rewriting Eq. (fTTj) as 



+ 3H(l + wf)p g = 0, 



we can define the effective equation of state for massive graviton 
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whose limit of a — > is — 1. When introducing the massive graviton £ maS s, both matter and 



massive graviton do not satisfy their own conservation law as is shown by ffT5]) and (p3 
Instead, the total mixture of matter and massive graviton fluid satisfies the conservation law 
as 

p t + 3H(p t +p t )=0, (20) 

with p t = p + p g and_»f_= p + p g . This shows clearly that our picture is quite different from 
([1]) imposed by Ref. 2j, yj. In other words, in order to express the evolution of universe due 
to the massive graviton properly, we will use its own native EoS Wg at and effective EoS w^ ff 
instead of u>de- This is a well-accepted approach to cosmology when two different matters 
coexist in the universe, showing that either decaying of massive graviton to CDM for Q < 
or decaying of CDM to massive graviton for Q > as in quintessence A(t)CDM model |{|, 
and the Brans-Dicke cosmology (jj]. 

Now we are in a position to solve Eq. <Q for a(t). For this purpose, Eq. ([9]) can be 
written as 



Hi 



+ + f (1-4 (2D 

where Hq is the Hubble constant at current, x is the current ratio of radiation density to 
dark matter density, and Q° m is the current density parameter of CDM. a g is the ratio of 
graviton mass and Hubble mass 



m r , 
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3.8026 x 10~ 6 %kg, 



(22) 
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where h Q is defined as H = 100h km/s/Mpc. Evaluating Eq. f l2T|) at the present time, the 
current scale factor is determined to be 



a 



In order to have real value for ao, one requires 



which corresponds to the inequality on the mass of massive graviton as 



m g > 2 v /l-n« m -xO> fl . 



Taking Q 



0.27 and x ~ 3.1 x 10 4 , we have two bounds on a g and m g 
a g > 1.71, m g > l.Hm H - 6.50 x 10- 69 /i kg. 



(23) 
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FIG. 1: a g as a function of current matter density parameter, Q^. Horizontal line represents its 
current value a g (fi^J = 1.71 at = 0.27. 



a g is depicted as a function of in Fig. 1 and its current value is a g (a^J = 1.71. 



At this stage, we introduce density parameters 
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Solving the Friedmann equations (121 j) together with the conservation laws (I15p and (I17p . 
the relevant cosmological parameters are determined to be 
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with a new variable x = ln(a/ao) G [—00,00] including 1 = at a = ao, instead of scale 
factor a. The time evolution of all parameters including density parameters, EoS, and 
Hubble parameter is shown as a function of x in Fig. 2. This graphs indicates that there is 
no phantom phase when using the native EoS u^ at and the effective EoS w e g but there is a 
phantom phase when using the geometric EoS u>de- This shows clearly that the description 
with u>de is not appropriate for interpreting the modified evolution due to the massive 
graviton. 




FIG. 2: Cosmological evolution for there density parameters [£l g (green) ,n m (magenta), f2 r (brown)], 
four equations of states [vjg at (blue) , u^ ff (cyan), u;DE(red), w eS (yellow 2j])], and Hubble parameter 
(black) as function of x, for = 0.27, x = 3.1 x 10 , and a g = 1.8. The bottom line corresponds 
to asymptote of —5/3 to u>de shown in Ref.^j, while the second lowest line represents asymptote 
of —1 to Wg S . 



In order to compute the age of the universe, we simply evaluate the integration 

r ao da 
'0 aH(a) 



(34) 



Substituting (1211 in fl3~4"|) leads to the age of the universe as shown Fig. 3 for chosen values of 
showing that the age of the universe is closely related to the mass of massive graviton. 
We note that H" 1 = lO.Ofy^Gy. 
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FIG. 3: Relation between the age tjj of the universe and the mass ratio a g for different values of f^, 
0.20(red), 0.25(blue), 0.27(green), 0.37(cyan), and horizontal line denotes 0.8, which corresponds 
to 8h 1 Gy. 

At this stage, we would like to mention that as is shown in Eq. (12T1) . the Hubble parameter 
becomes zero at a certain value of a c satisfying the condition of 



If one chooses a g as the saturating bound in Eq. fT24"|) . ao becomes zero from Eq. ff23l) and 
consequently, (1351) implies that a c becomes 1. However, this case is meaningless because the 
current scale factor is chosen to be zero. Hence a g should satisfy the bound of Eq. ( |24l) . 
This implies that the mass of massive graviton should be greater than a certain value to have 
a non-zero scale factor (m g > 1.71m#) at current time. In this case, we have the condition 
of a c > 1 which means that the energy density of massive graviton is negative as is shown in 
(|TT|) . In connection to this point, we may solve Eq. ( )2T|) for a as a function of cosmological 
time. The result is depicted in Fig. 4 for a g = 1.8, = 0.27, and x — 3.1 x 1CT 4 , indicating 
that there is no sizable difference when comparing to other cases. 

Finally, we derive the geometric equation of state wbe in ([I]) from our approach. We 




(35) 



S 



0.8- 




/ 0.6- 




0.4- 




0.2- 

i i i i i 





-2 



-1 



1 2 
t/t_H 



FIG. 4: Time evolution for scale factor a for a„ = 1.8, £1^ = 0.27, and x = 3.1 x 10 . 



start with the modified Hubble equation 



H 2 = Hi 
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Using x = ln(ao/a), we can rewrite the above two equations as 

8ttG 



where H 2 includes the contribution from all matter except CDM. Hence we can define 
'geometric' equations of state for CDM and other matters including massive graviton, re- 
spectively, through 

(37) 
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Differentiating Eq. (140]) with respect to x, we obtain 
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Rearranging this equation, one arrives at 
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which confirms that wde differs from the native EoS Wg at and equals to effective EoS Wg S 
when we use Eq. (j2T]l as Hubble parameter. Note that our wbe is different from that (2.14) 
of Ref. [3l even their forms are the same because their Hubble equations are different. In 

n 

addition, we can obtain effective EoS for total cosmological fluid as defined in Ref [2J: 

.3 . o? 2 
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which is depicted in Fig. 2. This EoS w eS is exactly the same as 

Pt 

w t = — 

Pt 



(43) 



(44) 



which is defined in (120)) . 

4. INTERACTION MECHANISM 

We have revisited the Massive Graviton Dark Matter scenario (MGCDM) which was 
recently considered as an alternative to dark energy models. When introducing the native 
and effective equations of state, it was shown that there is no phantom phase in the evolution 
toward the far past of a — > but the past accelerating phase appears. 




FIG. 5: Time evolution for interaction Q as a function of scale factor a with 1 its coefficient. We 
note that Q = at a = 0.74. 

How do we understand the past accelerating phase in the MGCDM? It is well known that 
the interaction between two different matters give rises to acceleration 10Nl2] . 

In the MGCDM, the interaction between massive graviton and CDM generates accelera- 
tion in the past. Graph of Fig. 5 shows that the CDM(p m , Q m ) decays to massive graviton 
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(pg, fig) for Q > and < a < 0.74, while massive graviton decays to CDM for Q < 
and 0.74 < a < 1. This picture is consistent with the behavior of (Q g ,Q m ) in Fig. 2. We 
note that the effective EoS Wg S is an decreasing function of x in the backward direction as is 
shown in Fig. 2. On the other hand, it is noted that the native EoS u>^ at is meaningless be- 
cause we are in the interacting picture. The graph of Fig. 6 shows the same graph depicted 
in terms of variable x. It shows that the CDM decays to massive graviton for Q > and 
— oo < x < 0.98, while massive graviton (p g , fl g ) decays to CDM (p m , fl m ) for Q < and 



0.98 < x < 1.16. The cosmic anti-fraction arisen 
graviton induces acceleration in the past universe 



rom the CDM decay process to massive 
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FIG. 6: Time evolution for interaction Q as a function of x with 1 its coefficient for = 0.27, 
X = 3.1 x 10 -4 , and a g = 1.8. We note that Q = at x = 0.98 and 2.95, and x = 1.16(o = 1). 



Finally, we would like to mention that the condition of positive energy density and nega- 
tive pressure ( fl3l) is not mandatory to be fixed. In Appendix, we can extend this condition 
by introducing the background static scaling factor a s . 

Consequently, any phantom phase des not appear in the evolution toward the far past. 
The interaction between massive graviton and CDM generates the past acceleration in the 
MGCDM. 
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Appendix: Background static metric with scaling a s 

If we choose the flat background metric go as 

ds 2 = c 2 dt 2 - a 2 [dr 2 + r 2 dnl] , (45) 

with a s an arbitrary scale factor, the relevant quantities are changed as 

3m 2 c 6 , 2 2 
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a 2 -a 2 ), (46) 
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We note that the positive energy density and negative pressure is required to explain the 
accelerating universe. In this case, one has the inequality 

a 2 < a 2 s . (48) 

Although the native EoS remains unchanged , the effective EoS w* s is changed as 



5a_ 

<(«) = — — (49) 



cff/ \ _ " 3? 



which is the same form as Eq. ( f!9l) except scaling factor a 2 . Also, the current scale factor 
is redefined by 



-^(i-K-xK)- (so) 

The condition of a g is slightly changed as 



<v B > ly/l-QP m - x Q° m , (51) 



a, 

which corresponds to the inequality condition on the mass of massive graviton as 

2 



m g >-y/l-n m -xn° m m H . (52) 
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That is, by introducing the scaling a s of the background static metric, it is possible to restrict 
the maximum scale factor for the universe as in fj48l) . 
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